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Abstract — A Lyapunov-based control design for natural trajectory- 
tracking problems is analyzed for quantum states where the analysis 
in the generic case is not applicable. Using dynamical systems tools we 
show almost global asymptotic stability for stationary target states subject 
to certain conditions on the Hamiltonians, and discuss effectiveness of 
the design when these conditions are not satisfied. For pseudo-pure target 
states the effectiveness of the design is studied further for both stationary 
and non-stationary states using alternative tools. 

I. Introduction 

Recent technological advances have prompted significant interest 
in developing the foundations of quantum control theory. One major 
concern is how to design the dynamics to steer the system state 
to a stationary target state, or track the natural trajectory of non- 
stationary states. Numerous results have been obtained on the theory 
of Lyapunov-based control [1]-[11] and the references in [11]. Most 
papers have used the Hilbert- Schmidt (HS) distance as the most 
natural Lyapunov function and in this setting [4]-[6] showed that 
the target state is asymptotically stable under a sufficient condition 
equivalent to controllability of the linearized system for pure states 
represented by wavefunctions, although an additional control had to 
be added to fix the relative phase of the state. An alternative design 
based on an implicit Lyapunov function was proposed in [7] to render 
(pure) target states asymptotically stable when the linearized system is 
not controllable. The more general case of systems whose states must 
be represented by density operators was recently considered in [9], 
[10], but our analysis for generic quantum states [12], for instance, 
showed the dynamical landscape and convergence behavior to be 
more complicated than described in [10]. The aim of this TechNote 
is to extend this analysis to non-generic states, which, although a 
set of measure zero, deserve special consideration as they include 
important cases such as pseudo-pure states. 

Specifically we show that, although the LaSalle invariant set for 
non-generic states is much larger than in the generic case, for ideal 
systems the target state not only remains isolated and thus locally 
asymptotically stable, but a dynamical systems analysis for stationary 
target states shows that all other critical manifolds are unstable and 
attract at most a measure-zero subset of the state space. At least for 
the stationary target states, these arguments rigorously show that we 
indeed have almost global convergence as claimed in [10], although 
the set of states not attracted to the target state is significantly larger 
than claimed in [10]. The conditions for effectiveness of the method 
for stationary target states are similar to those for generic states, and 
we again find that center manifolds emerge around the target state 
when these conditions are not satisfied, but in contrast to the generic 
case, the emergence of a center manifold now is conditional on the 
target state, i.e., the method may still be effective for some target 
states even if the ideal Hamiltonian criteria are not satisfied. For the 
special class of pseudo-pure states, we can further show explicitly that 
for systems with ideal Hamiltonians, the method is effective not only 
for stationary target states but also for almost all non-stationary ones. 
This is similar to the situation for the generic case but in this special 
case we can characterize the non-regular target states explicitly, which 
shows that they include some states of special interest such as so- 
called CAT states and maximally entangled states. 



II. Lyapunov control and LaSalle invariant set 

We study the bilinear Hamiltonian control problem for a quantum 
system on an n-dimensional Hilbert space 7i. The state of such a 
system is generally represented by a trace-one positive operator p on 
TL, referred to as density operator, whose evolution is governed by 
the quantum Liouville equation 

p(t) = -i[H + f(t)H 1 ,p(t)], (1) 

where f(t) is an admissible real- valued control, and Ho and Hi are 
time-independent Hamiltonians corresponding to free evolution and 
control interaction terms, respectively. 

Definition 11.1. A density operator p represents a pure state if it is 
a rank-one projector, and a mixed state otherwise; p is generic if it 
has n distinct eigenvalues, and pseudo-pure if its spectrum has only 
two distinct eigenvalues occurring with multiplicities 1 and n — 1. 

The control problem is to design a control f(t) such that the system 
state p(t) converges to the target state pa as t — > oo. For Hamiltonian 
evolution a necessary condition for pd to be reachable from an initial 
state p(0) is that both p(0) and pd have the same spectrum, and we 
shall assume this condition to be satisfied here. The set of density 
operators isospectral with pd forms a compact manifold A4, whose 
dimension depends on the spectrum of pd- Pd is a stationary state of 
the system if and only if [Ho, pa] = 0. If pd is not a stationary state 
then the control problem becomes a trajectory-tracking problem, and 
to be able to apply LaSalle's invariance principle to this case, we 
formulate the control dynamics on an extended state space M x M: 

pit) = -i[Ho + f( P , Pd )H u p{t)l (2a) 
p d (t) = -i[H ,p d (t)], (2b) 
f(p,p d ) = Tr([-iHi,p(t)]p d (t)), (2c) 

where the control function / is chosen such that 

V(p, p d ) = | Tr[(p - pd) 2 } = Tr(p d ) - Tr(pp d ) (3) 

is non-increasing along any flow, i.e., 

f t V(p, Pd ) = -f(t)Tr(pd[-iHi,p}) = -f(t) < 0. 

Thus, the dynamical system (p(t), pd(t)) G M x M is autonomous, 
all solutions are bounded, and V(pi,p2) is a Lyapunov function on 
A4 x M. From LaSalle's invariance principle [13] we have: 

Theorem 11.1. Any evolution (p(t), pd(t)) for the dynamical system 
|2]> will converge to the LaSalle invariant set, E = {(p\,p2) G 
M x M\V(p(t), pd(t)) = 0, (p(0), Pd (0)) = (pi,p 2 )}. 

The effectiveness of the control in steering the system towards 
the target state depends on the asymptotic stability of pd- Local 
asymptotic stability is a prerequisite for almost global convergence, 
and a necessary and sufficient condition for the former is that the 
target state be isolated in E. Hence, we first characterize the invariant 
set E, which depends on both the Hamiltonians Ho, Hi and p d . 

Definition II.2. Choose a basis such that Hq = diag(ai, . . . , a n ), 
which is always possible as Ho is Hermitian. Let Hi = (b^i) and 
Uki = at — ai be the transition frequencies. System ^ is ideal 

(i) Ho is strongly regular, i.e., ui^i 7^ u pq unless {k,£) — (p,q). 

(ii) Hi is fully connected, i.e., bte 7^ except (possibly) for k = £. 

As shown in [12], we have 

Theorem II.2. The invariant set of ([2]) for an ideal system satisfies 
E = {(pi,P2) G M x M : [pi,p 2 ] = diag(ci, . . . ,c„)}. 
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III. Stationary target state 

For a generic stationary state pd, it has been shown that the LaSalle 
invariant set E is comprised of n\ distinct stationary states, which 
coincide with the n! critical points of V(p) = V(p,pd) as a Morse 
function [12]. If pd is not generic then E is much larger and the 
topology becomes quite complicated. If pd is stationary then j2| can 
be reduced to an autonomous system on M, and the LaSalle invariant 
set to E — {p : [p, pd] ~ 0}. We can also choose a basis such that 
Ho and pd are simultaneously diagonal, Ho — diag(ai, . . . , a n ), 



p d = diag(«;i, . 



(4) 



where {we }e=i are the distinct eigenvalues of pd with multiplicities 
ne and ne = n. The state space in this case is the flag manifold 

M ~U(n)/[U(ni) x ... xU(n fc )] (5) 

of dimension iim = n' 2 — X/f-i n f ■ 

Definition III.l. xo is a critical point of f(x) if^f(x)\ x = xo = 0; 
it is hyperbolic, if the Hessian at the critical point is nonsingular. f 
is a Morse function if all of the its critical points are hyperbolic. 

For a given pd, we can investigate the critical points of the 
Lyapunov function V(p) = V(p,pd) on M. 

Lemma III.l (proof in [12]). For a given stationary target state pd 
the critical points ofV(p) = V(p, pd) on M are such that [p, pd] = 
0. 

If pd is generic, there are n! critical points of V in total; if pd is 
non-generic, there are more critical points, forming different critical 
manifolds (see also [14]). 

Example III.l. For a three-level system with pd = diag(|, \,\) 
the set of p G M — {pd} that commute with pd has the form 



po = 



I Hill U>12 N 

W*2 l«22 
V i, 



which forms a manifold Mo isomorphic to the Bloch sphere, on 
which V assumes its global maximum, while p — pd corresponds to 
minimum of V . 

Theorem III.l. Let pd be a given stationary target state. Then p = 
Pd is an isolated hyperbolic critical point of V(p) corresponding to 
the global minimum and there are no other local or global minima. 

Proof: V(p) clearly assumes its global minimum for p — pd- 
To show p is a hyperbolic minimum of V it suffices to show that 
it is a hyperbolic maximum of J(p) = Tr(ppd). Choose a basis so 
that pd is diagonal and let p — po be a critical point. Any point 



in a neighborhood of po can be written as p 



poe 



where 



3 — {Xke, Xke, Afe} is the basis of the Lie algebra su(n) defined in 
the Appendix. Substituting this into J gives J = Tv(e x " poe~ x a Pd)- 
To show that the critical point po = pd is a maximum of J, we need 
to find Tim — dim.M independent directions in which J is a local 
maximum. If we choose curves through po with x • a = X k et then 

J = Tr(poPd) + i 2 {- Tr(p A^p d A fc£ ) 

+ | Tr(poPdXtt) + | Tr(p () ALpd)} + e(|i| 3 ). 

The conjugate action of Xke on p = pd on the (k, £) subspace swaps 
the fc-th and £-th diagonal elements. If we choose the curve with 
x • a = Xke, we get a similar result. Hence the number of swaps that 
decrease the value of J is 
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n e -\ + n t _in t 



k 

E 

<=i 
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dim M. 



Thus p = pd is a hyperbolic point of J, which is necessarily isolated. 
At other critical points, it is easy to see that there always exists some 
swap Xke that increases the value of J, and hence all the other critical 
points cannot be local maxima of J, or local minima of V . ■ 

Remark III.l. For a given pd, the critical points ofV(p) = V(p, pd) 
form a finite number of isolated critical manifolds with pd as an 
isolated hyperbolic minimum. When pd becomes generic the critical 
manifolds are reduced to n\ isolated hyperbolic critical points. 



From Lemma |III.1| we can see that for ideal Hamiltonian and 
stationary p<j, the LaSalle invariant set E coincides with the critical 
points of V(p) = V(p, pd), and there are p — ni " l nk] stationary 
solutions including pd, either isolated or located on an isolated 
manifold in E. In particular, since pd is isolated in E, the LaSalle 
invariance principle guarantees its local asymptotic stability, but if 
we want to investigate the asymptotic stability and convergence 
properties on a large scale, we must analyze the eigenvalues of the 
linearized dynamics at the other stationary states. 

Theorem III. 2. For ideal systems the stationary state pd is a hyper- 
bolic sink of the dynamical system ([2j and thus locally asymptotically 
stable, while all other stationary solutions have unstable manifolds. 
Hence, pd is almost globally asymptotically stable. 

Proof: In order to analyze the stability at the stationary state we 
reformulate the dynamical system ([2j in the Bloch representation 

k(t) = (Ao + f(s, s d )A!)s(t) (6a) 

k d {t) = Aos d (t) (6b) 

f(s,s d ) = Sd T Ais, (6c) 
2 —i 

where s — (sfc)£ =1 with = Tr(afep), and Ao and Ai are the 
anti-symmetric matrices 



Aa(k,k') = Tr(iHo[a k ,a k ,}), 
A 1 (k,k') = Tr(iH 1 [a k ,a k ,]), 



(7a) 
(7b) 



<Jk being the elements of the basis of su(n) denned above. For 
stationary pd, the dynamics is reduced to 



(8a) 
(8b) 

(9) 



s(t) = (Ao + f&Ajm 

f(s) = 3d AlS. 

The linearized system near the critical point so is 

s = Df(s ) ■ (s- s ), 

where Df(so) — Ao + AiSo ■ s~d T Ai is a linear map defined 
on K n _1 . (The map is denned on R n _1 as we neglect the so 
coordinate of p, which is constant due to trace preservation.) In 
order to show s d is hyperbolic, it suffices to show that there are 
tim eigenvalues with nonzero real parts, corresponding to um 
eigenvectors in the tangent space of M at Sd, denoted as TM(sd). 

2_i 

Let Sc and St be the subsets of R™ corresponding to the Cartan 
and non-Cartan subspaces C and T of su(n) = C ® T. Since 
[— iHo,pd] = and [— iHi,pd] G VT , we have AoSd = and 
j4iSd G St- Let Hi = (b k e), and v be a column vector consisting 
of |n(n — 1) blocks 

-4k,e) _ Afc^ / 9(6fef ) 
~ V2 \®(bkt) 

Let B = Bo — vv T be the restriction of Df(sd) to the subspace St 
as before. Following a similar argument as in [12] it is easy to see 
that for (k,£) such that = 0, the eigen-element (uike, &ke) of Bo 
is also an eigen-element of B as v T e k e = 0, and that e k e corresponds 
to a direction orthogonal to the tangent space Tj^\(sd). The number 



(10) 
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of such (k,£) is TV = 2j2t=i ("'). We can therefore show that 
the remaining eigenvalues of B with eigenvectors corresponding to 
the directions in TVi(.Sd) must have non-zero real parts. A simple 
counting argument shows that the number of these eigenvalues is 
2(2) — n = dim(.M) and thus pd is a hyperbolic point. Since pd 
achieves the minimum of V, these eigenvalues must have negative 
real parts, i.e., pd must be a sink. Hence, there exists a neighborhood 
M of pd such that p(0) £ N will converge to pd for t — > +00, 
which establishes local asymptotic stability of pd. 

We can similarly show that any other stationary state po that is 
isolated in E is a hyperbolic fixed point, i.e., all eigenvalues of 
linearized system at po have non-zero real parts. As po is also an 
unstable critical point of the Lyapunov function V by Theorem |III.1| 
we can conclude there is an unstable manifold at po, as in the 
generic case [12]. When the stationary state po lies on the critical 
manifold in E then the linearized system at po must again have an 
unstable manifold; for otherwise all eigenvalues of Dj at po would 
all have non-positive real parts, and p would be a local minimum 
of V, contradicting Theorem |III. 1 1 Therefore, all stationary states 
except pd have an unstable manifold of positive dimension, and the 
dimensions of all stable manifolds at these unstable stationary states 
are less dim.M. Hence, almost all non-stationary solutions will not 
converge to any of these unstable stationary states or their center 
manifolds, and therefore must converge to pd- Thus, pd is almost 
globally asymptotically stable. ■ 

Example 1 1 1.2. With pd = diag(j, |, «) as in the previous example 
we have um = 3 2 — 2 2 — 1 = 4 and E contains pd and two other 
stationary states pi = diag Q, |, and p2 = diag (|, |, |). 
Analysis of the linearized dynamics shows that the two tangent 
vectors of the center manifold at pi are also the tangent vectors 
of E. Therefore, except for pd, which is isolated, the points in E 
form the center manifolds at the stationary states pi, i — 1,2. 

For generic target states we showed in [12] that for non-ideal 
systems the target state itself generally becomes a center on an 
attractive center manifold, and the control becomes ineffective. For 
example, when Ho is regular but not strongly regular then the 
LaSalle invariant set E usually becomes much larger, forming a 
center manifold around the target state pd, and most solutions p(t) 
converge to points on this center manifold other than pd- This is 
still true for most non-generic stationary pd, for the same reasons, 
but unlike in the generic case, there are special target states that are 
asymptotically stable even if the Hamiltonian does not satisfy the 
criteria for ideal systems. 

Example III.3. Consider the three-level system with Hi fully con- 
nected hut Ho not strongly regular 



Ho = 



Let pd = diag(ai, . . . , a n ) be the stationary target state, p — (/3 mn ) 
and Aki = Qfc — a.£. Then 

/ -/?i 2 Ai2 -ftsAasX 

[p, Pd )= [jSJjAia -ftsAas , 

\fit 3 A 2 3 /%A 2 3 / 

but for the given Hamiltonians Ho and Hi, we find that any point p 
in the LaSalle invariant set satisfies 
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Fig. 1: Evolution of V(t) — V(p(t), Pd(f)) with j/-axis in logarithmic 
scale. Each graph shows V(t) for N = 50 different initial states p(0) 
for a model system with Hi fully connected and Ho regular but not 
strongly regular, (a) For pd = diag(l,0,0) all trajectories appear 
linear with negative slope in the semi-logarithmic plot, indicating 
exponential convergence; (b) for pd = diag(0, 1, 0) most trajectories 
flat-line after some time at different non-zero distances, suggesting 
those trajectories do not converge to pd but its center manifold. 



If pd — diag(0, 1,0), then any p G E must satisfy fin = 
fivx = 023, |/3ia| = fill and \fii 2 \ = fin - ^fiii- Analogous to 
the generic case, the points in LaSalle invariant set form a center 
manifold around pd and any solution starting outside E will stop 
converging to pd after certain time, as shown in Fig. [7] f b ). 

However, if pd = diag(l,0,0) then any p G E — {pd} satisfies 

/0 \ 

p = fi 2 2 P23 
\0 fih fiuj 

which corresponds to the maximum value ofV(p, pd). Therefore, any 
solution starting outside E satisfies V(p(0), pd) < Knax and hence 
converges pd, as shown in Fig. \7Ja). 

This example shows that the non-generic case differs from the 
generic one, but for most non-generic target states, a center manifold 
around the target state will still appear for non-ideal Hamiltonians, 
rendering the control design ineffective. 

IV. Pseudo-pure target States 

Pseudo-pure states, although a subset of measure zero of all 
possible states, deserve special consideration as they include the 
important special class of pure quantum states, and non-pure pseudo- 
pure states play an important role in applications such as ensemble- 
based quantum information processing. Moreover, for pseudo-pure 
target states the analysis simplifies, allowing us to derive the strong 
stability results for both stationary and non-stationary target states. 

Any pseudo-pure state pd(0) can be written 



p d (o) = W n + 1 — f n ( f = wH + u n,| 

71—1 



<w < 1, (11) 
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where Ho = |^ r o){^ r o| is a rank-1 projector onto some pure state 
j*&o), and is the projector onto the orthogonal subspace. 

Theorem IV. 1. Assuming ideal Hamiltonian, any solution p(t) such 
that V (p(0) , pd(0)) < Vmax converges to Pd(t) as t —> +00 except 
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when pd(0) = (r k e) has a single pair of non-zero off-diagonal entries 
of the form r k e = |(w — u)e' a and r kk = ru = §(n> + «)• /» f^e 
exceptional case any solution p(t) converges to the orbit of Pd(t) but 
in general p(t) -/-> pd(t) as t — > +00 anrf V(p,pd) can take any 
limiting value between and Vm ax . 

Proof: For any (pi,p2) £ -E, pi and p2 must also be pseudo- 
pure, with the same spectrum {w, u}, i.e., p k = wll k + ull k for 
k = 1,2, and we have 

[pi,Pa] 



( m - u ) 2 [ni,n 2 ]. 



(12) 



Thus the LaSalle invariant set contains all points such that M = 
[IIi. II2] is diagonal, according to Theorem II.2 Let U k = |^ , ^)(^ r fc|, 
k = 0, 1, 2. Setting 

= (die \...,a„e ") , 
1*2) = (61 e M1 , 



, b„e ' 



we have M = |*i)(*i|*2}(*2| - *2}(*2|*i)(*i|. For 
{pi,p2) £ E, the off-diagonal components of M must vanish: 

Mm = a fc 6ie i(a *- ft) (* 1 |* 2 >-a^ fc e i(3 '=- a * ) (* 2 |*i) = 0, (14) 

for all 7^ £. Let (*i|* 2 ) = re lS . We have the following two cases. 

(a) r = i.e (*i|* 2 ) = 0. In this case, [pi,p 2 ] = 0, and 

V(pi, p 2 ) = Vmax = (W - U) 2 . 

(b) If r then together with M fcfc = gives: 

.a k b e — a e b k , ft + f3 e — a k + a e + 28 (15) 

If a k — then = a k bg = atb k for £ 7^ k and we must have b k = 
as a,i = is not allowed as a is a unit vector. Ditto for b k — 0. 
Let 7+ be the set of all indices k so that a k , b k 7^ 0. Then 



— = — , Vk,£el+ 
b k be 



(16) 



and thus a — yb, where a = {ak) and b = (b k )- Since a and b are 
unit vectors in R" , 7 = 1 and a — b. 

As for the phase equations l |15[ >, if = fefc = then Mki = is 
automatically satisfied, thus the only non-trivial equations are those 
for k,£ £ /+. If the set J+ contains m > 2 indices then taking 
pairwise differences of the m(ni — l)/2 non-trivial phase equations 
and fixing the global phase of \^k) by setting a ni = f3 ni — shows 
that a — f3. For example, suppose I + — {1,2,3} then we have 3 
non-trivial phase equations 

ft + (5 2 = qi + a 2 + 26», 

Pi + = ai + a 3 + 26», 

ft + #3 = «2 + Q 3 + 26>, 

taking pairwise differences gives 



ft 


-ft 


= «2 


— «3 


ft 


-ft 


= ai 


— C*3 


ft 


-ft 


= Oil 


— Q?2 



and setting 03 = ft = shows that we must have 02 = ft and 
0.3 — ft. Thus, together with a — b we have pi = p 2 . If /+ 
contains only a single element then and |$ 2 } differ at most by 
a global phase and again pi = p 2 follows. Incidentally, note that for 
|*i) = |* 2 ) we have (*ii* 2 ) = 1, i.e., r = 1, 9 = 0. 

The only exceptional case arises when 7+ contains exactly two 
elements, say {1,2}, as in this case there is only a single phase 
equation ft + ft — ai + a.2 + 28, and thus even fixing the global 
phase by setting 012 — ft = 0, only yields ft — ai = 28. This 
combined with a — b gives 

iS ;, T , I.T, \ 2 2i6 , 2 

re = (Vi|W 2 ) = a 1 e + a 2 




Fig. 2: V^t) = V(p(t), Pd(t)) with y-axis in logarithmic scale. Each 
graph shows V(t) for iV — 50 different initial states p(0) for a two- 
qubit system with ideal Hamiltonian. (a) For the target state |$) = 
^g(l,0,0,2) T all trajectories converge exponentially to the target. 



(b) For the target state | *) 



(1,0,0, 1) , most trajectories stop 



converging after some finite time and lim^oo V(t) > 



and thus a\e iB + a\e~ %e = r or 2i sin 8(a'( - a|) = 0. Therefore, 
either 8 = or ai = a 2 . If 6* = then a = /3 and pi = p 2 , which 
is one possible solution in E. If 8 7^ 0, then any (pi,p 2 ) satisfying 



o) J 



|*i) = 2- 1/2 (l,e"\0, 
!* 2 ) = 2- 1/2 (l,e l,3 ,0,...,0) 



(17a) 
(17b) 



with j3 — a — 28 is also in E. Hence, if |^o} has only two nonzero 
components with equal norm, e.g., if 



Mo) = 



Vo 



ri 2 
rn 











u 







(18) 



with ru = \{w + it), ri 2 = \(w - u)e la , and |* ) = 
2 1/2 (1, e ia ,0, . . . , 0) T , then E contains all points (pi,p 2 ) satisfy- 
ing l |17[ l, which includes pi = p 2 and pi _L P2. Since pi and P2 
lie on the orbit of pd(0), any solution p(t) will converge to this 
orbit but we cannot guarantee p(t) — + Pd(t) as f — > +00. For all 
other pd(0) E contains only points with either pi = p2 or pi _L p2, 
corresponding to V = and V = Vmax, respectively, and since V is 
non-increasing, any solution p(t) with V(p(0), p<j(0)) < Vmax will 
converge to pd(t) as t — > +00. ■ 

V. Applications for two-qubit systems 

Four-level (or two-qubit) systems are interesting for many reasons, 
but for our purposes four-level systems are of interest as n = 4 is 
the lowest dimension that admits states that are neither generic nor 
pseudo-pure. Any stationary state pd has one of the following forms: 

• generic: pd = diag(a, b, c, d) 

• pseudo-pure: pd = diag(a, b, b, b) 

• other: pd = diag(a, a, 6, b) 

where a, b, c, d are distinct. The first two cases correspond to the cases 
discussed above and in [12]. In the third case, pd = diag(a, a, b, b), 
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the state manifold M is homeomorphic to the flag manifold 
U(4)/U(2) x U(2), or SU(4)/{exp(a) : a g C®T U ®% A } and has 
real dimension tim — 8. Let (a, a, b, b) represent diag(a, a, b, b), etc. 
There are six stationary solutions: (a, a, b, b), (a, b, a, b), (a, b, b, a), 
(b,a,a,b), (b,a,b,a) and (b,b,a,a). According to the results in 
Section III the critical points of V(p) = V(p,pd) consist of 
six components: two isolated hyperbolic critical points, (a, a, b, b) 
and (b,b,a,a), corresponding to the global minimum V = and 
maximum V = 2 (a — b) 2 , respectively, and four critical manifolds 
homeomorphic to S 2 x S 2 on which V assumes the intermediate 
value (a — b) 2 . 

Calculating the Df(so) restricted to St from the linearized 
dynamics l|9j at each stationary state po, as in the proof of Theo- 
rem |IIL2] shows that for ideal systems, the dimension of the stable 
manifold at pd = diag(a, a, b, b) is n s — 8 = tim, while at 
(6, b, a, a), the dimension the unstable manifold n u — 8 = tim- Thus 
these two points are a hyperbolic sink and source, respectively, with 
Knin = and Vmax = 2(a — b) 2 . At any of the other four stationary 
states, i.e., (o, b, a, b), (a, b, b, a), (b, a, b, a) and (b, b, a, a), we have 
Vo = {a — b) 2 , a stable and unstable manifold of dimension 2 
and two pairs of purely imaginary eigenvalues, showing that these 
intermediate critical points are centers, but the dimension of the 
corresponding center manifolds is less than um, and each of these 
centers has an unstable manifold. Therefore, the center manifolds and 
the set of states that are attracted to them, form a set of measure zero 
in the state space, and even in the vicinity of a center manifold, many 
states will converge to pd- Thus we can say that pd is almost globally 
asymptotically stable. 

While this suggests that the Lyapunov design is mostly effective, 
at least for ideal Hamiltonian, unfortunately, it still fails for some 
of the most interesting problems. As an example, consider the 
problem of Bell state generation, i.e., tracking the trajectory of a 
maximally entangled state such as |^+) = -^j(|00) + 1 11)), which 
is important in quantum information. Setting pd(0) ~ \^+){9+\ 
we see immediately that we cannot track this state because pd(0) 
satisfies the conditions of Theorem IIV. 1 1 for which the method fails: 
Since pd(0) is pure we have w = 1 and u — and taking 
the standard Pauli product basis we have pd{0) = {vkl) with 
rn = ri4 = V4i — r-44 — ~ and r k i = otherwise. Fig.pfb) shows 
that most initial states converge to points in the invariant set E that are 
finite distance from pd- We could console ourselves that the method 
does work well for most other target states such as the example in 
Fig. |2ja). However, the effectiveness of the method in these cases is 
predicated on the assumption of ideal Hamiltonians, which excludes 
most typical physical Hamiltonians such as H — 0.1 er z (g> a z and 
Hi = a x ® I + 0.9 1 (g> a x for a system with constant Ising coupling 
and local 2>rotations. In this case Ho is not even regular and many 
of the off-diagonal elements of Hi are zero, and the invariant set E 
becomes huge, rendering the method ineffective. 

VI. Conclusion 

We have extended our previous work on Lyapunov control for 
generic quantum states to the non-generic case. Although the LaSalle 
invariant set is now much larger, containing critical manifolds, we find 
that for ideal Hamiltonian, the stationary target state is still almost 
globally asymptotically stable. For non-ideal systems we find that 
in general the target state changes from an isolated asymptotically 
stable minimum to a center on a center manifold as in the generic 
case, but unlike in the generic case, there are special cases where a 
non-generic target state remains asymptotically stable under non-ideal 
Hamiltonian, though these states are exceptional. For the important 
class of pseudo-pure states, necessary and sufficient conditions for 



convergence to the target state for both stationary and non-stationary 
cases are derived by alternative mean. Finally, application to two- 
qubit systems shows that even for ideal systems the small set of 
(non-stationary) states that are non-trackable using Lyapunov control 
includes some of the most interesting cases such as CAT states and 
maximally-entangled Bell states. 
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Appendix 

A standard basis for the Lie algebra su(n) is given by {X k e, X k e, X k } 
for 1 < k < t < n, where 

Afc = i(e kk — ek+i,k+i) 
Am = i{eu + eik), ^ki = {&kl — e«) 



and the (k,l) entry of the matrix e mn is 8kmSe n , and i = 
We have the useful identities 

Tr(\kt\ k >t>) = Tr(\ M \ kU >) = -25 k y&u> (20a) 
Tr(A M W) = (20b) 

and for a diagonal matrix D — ^2-1 dk&kk we have [D, A*,] = 0, 

[D, X ke ] = +i{d k - d e )X ke , [D, X ke ] = -i(d k - d e )X ke . (21) 

The basis {19} is not orfhonormal but we can define an equiv- 
alent orthonormal basis by normalizing the n 2 — n non-Cartan 
generators -^X k i and -^X k e, and defining the n — 1 orthonor- 
mal generators for the Cartan subalgebra o n i_ n+r = i[r(r + 
!)r i/2 {YZ=i - re r+ i, r+ i) for r = l,...,n- 1. 
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